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Abstract 

Utilization of "hidden genes" from wild species has emerged as a novel option for enrichment of genetic diversity for 
productivity traits. In rice we have generated more than 2000 lines having introgression from 'A' genome-donor wild species 
of rice in the genetic background of popular varieties PR1 14 and Pusa44 were developed. Out of these, based on agronomic 
acceptability, 318 lines were used for developing rice hybrids to assess the effect of introgressions in heterozygous state. 
These introgression lines and their recurrent parents, possessing fertility restoration ability for wild abortive (WA) cytoplasm, 
were crossed with cytoplasmic male sterile (CMS) line PMS17A to develop hybrids. Hybrids developed from recurrent 
parents were used as checks to compare the performance of 318 hybrids developed by hybridizing alien introgression lines 
with PMS1 7A. Seventeen hybrids expressed a significant increase in yield and its component traits over check hybrids. These 
17 hybrids were re-evaluated in large-size replicated plots. Of these, four hybrids, viz., ILH299, ILH326, ILH867 and ILH901, 
having introgressions from O. rufipogon and two hybrids (ILH921 and ILH951) having introgressions from O. nivara showed 
significant heterosis over parental introgression line, recurrent parents and check hybrids for grain yield-related traits. Alien 
introgressions were detected in the lines taken as male parents for developing six superior hybrids, using a set of 100 
polymorphic simple sequence repeat (SSR) markers. Percent introgression showed a range of 2.24 from in O. nivara to 7.66 
from O. rufipogon. The introgressed regions and their putative association with yield components in hybrids is reported and 
discussed. 
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Introduction 

Rice is grown annually on more than 153 million hectares with 
an annual production of about 483 million tons globally [1]. 
Remarkable progress has been achieved in rice production and 
productivity during the past 50 years, but still, food security is 
among the biggest challenges of the 21st century. Global rice 
demand is estimated to rise from 676 million tons in 2010 to 763 
million tons in 2020 and further increase to 852 million tons in 
2035. This is an overall increase of 26% or 176 million tons in 25 
years [2]. Rice productivity witnessed a tremendous increase 
during the green revolution phase (late 1960s to 1970s); however, 
these gains slowly reached a plateau. Hybrid-rice technology is 
advocated as a viable option for improving productivity, as has 
been amply demonstrated in Peoples' Republic of China [3]. 
Tremendous genetic variability exists in cultivated rice, but this is 
just five percent of the total variability existing in the genus Oryza 
[4] . Plant breeders have recognized the worth of wild species and 
used these for improving simply inherited traits. The most 
successful examples of utilizing wild species include the use of O. 
spontanea as a source of wild abortive cytoplasmic male sterility [5] , 
0. nivara genes to provide resistance against grassy stunt virus [6] , 



and 0. longutaminata gene Xa21 [7]- [8], 0. rufipogon gene Xa23 [9] 
and 0. nivara gene Xa38 [10] for resistance against bacterial blight. 
There is growing concern that limited genetic diversity in present 
day high-yielding rice germplasm impedes further improvement in 
productivity. Incorporation of potentially favorable alleles from 
wild ancestors of rice into improved genotypes for productivity 
traits has emerged as a promising strategy. The foundation of 
utilizing productivity-enhancing genes from wild species was laid 
by Frey et al. [11] in oats. Following the availability of advanced- 
backcross quantitative trait loci (AB-QTL) approach proposed by 
Tanksley and Nelson [12], discovery of favorable genes from wild 
species received impetus as a plant breeding strategy. Conse- 
quently, several yield-enhancing alleles from wild species of rice, 
such as 0. rufipogon, were identified and introduced into high- 
yielding elite cultivars [ 1 3] — [2 0] . Apart from 0. rufipogon, search for 
yield-enhancing quantitative trait loci (QTL) alleles was success- 
fully extended to 0. glaberrima [21], 0. grandiglumis [22], 0. 
glumaepatula [23], 0. longistaminata [24] and 0. minuta [25]. The lines 
generated provided precise estimates of genetic effects of 
introgression under relatively uniform and elite genetic back- 
grounds [12]. Furthermore, there are several indications that 
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Figure 1. Development of alien introgression lines. 

doi:10.1371/journal.pone.0096939.g001 



introgression lines served as dynamic entities to identify new genes 
[26], distinguish pleiotropy from linkage [27] and identify QTL 
contributing to heterosis, particularly those showing overdomi- 
nance [28]. 

In rice, superiority of introgression lines over recurrent parents 
has been reported by many workers, but the effect of alien 
introgression in F] hybrids has not been studied extensively. This 
study was conceived with the hypothesis that alien introgression(s) 
in homozygous condition may not be able to contribute positively 
to the phenotype because of replacement of a cultivated genome 
segment. However, in heterozygous condition, the recipient 
complement is conserved and minor negative effect associated 
with the alien segment is likely to be masked. Therefore, the 
present study was planned to determine the effect of novel genes 
from wild species on the yielding ability of hybrids and to locate 
the alien segments responsible for enhanced magnitude of 
heterosis in rice hybrids. 

Materials and Methods 

Development of introgression lines by backcross 
breeding 

A backcross breeding programme was initiated to introgress 
useful variability for improving yield and yield components from 
wild Oryza species into two high yielding indica rice varieties, viz., 
Pusa44 and PR 114. These varieties are widely adapted and restore 
fertility in wild abortive (WA) cytoplasm from O. spontanea. Forty- 
five accessions of various 'A' genome wild species, viz., O. rufipogon, 
0. nivara, 0. barthii, O. longistaminata, 0. glumaepatula, 0. meridionalis 
and cultivated African rice O. glaberrima were used as donors. Two 
to three backcrosses, followed by continuous selfing for five or 
more generations, led to the development of >2000 alien 
introgression lines (ILs). During backcrossing and selfing, pheno- 
typic selection was practiced for yield, yield component traits and 
phenotypic acceptability (Fig. 1). 



Development and evaluation of hybrids 

Based on agronomic evaluation of a large set of introgression 
lines in an earlier experiment, 318 ILs were selected and crossed 
with a CMS line PMS 1 7 A, carrying 'WA' cytoplasm to develop 
test hybrids. Among the 318 ILs used for generating test hybrids, 
99 were derived from 0. rufipogon, 99 from 0. nivara, 79 from 0. 
glaberrima, 30 from O. longistaminata, six from O. glumaepatula, five 
from 0. barthii and two from 0. meridionalis, which were used as 
donors. Recurrent parents, viz., Pusa44 and PR114 were also 
crossed with the same CMS line to develop check hybrids. Both 
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Figure 2. Development of Test and Check hybrids. 

doi:10.1371/journal.pone.0096939.g002 
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test and check hybrids had similar genetic background except for 
the alien introgression segments transferred from wild species 
(Fig. 2). 

In 2008, test hybrids, check hybrids and recurrent parents were 
evaluated at Punjab Agricultural University Ludhiana, India for 
yield and yield-contributing traits in a simple lattice design with 
two replications and plot size of 1.26 m 2 . Observations were 
recorded on days to 50% flowering, days to maturity, plant height 
(cm), percent pollen fertility, panicle length (cm), spikelets per 
panicle, grains per panicle, percent spikelet fertility, 1000-grain 
weight (g) and grain yield (t/ha) following the standard evaluation 
system adopted at IRRI [29]. Based on plot yield and other 
contributing traits, 1 7 promising hybrids were identified. F! seeds 
of the 17 promising test hybrids and two check hybrids, viz., 
PMS17A/Pusa44 and PMS17A/PR1 14, were generated at 
Central Rice Research Institute (CRRI), Cuttack, India during 
winter season 2008-2009 for further evaluation. In 2009, the 1 7 
test hybrids along with hybrid checks and recurrent parents were 
evaluated at Ludhiana in a randomized block design (RBD) using 
a plot size of 7.02 m 2 for each test entry. This evaluation was 
carried out in two experiments. Experiment I was planned to 
evaluate four test hybrids developed from Pusa44-derived ILs and 
13 test hybrids developed from PR1 14-derived ILs were evaluated 
in experiment II. Statistical analyses were performed on the basis 
of entry means in each replication. Analysis of variance for simple 
lattice and randomized block design was carried out using software 
SPAR 2.0 [30]. 

Marker analysis of introgression lines 

DNA from the ILs, recurrent and donor parent lines was 
isolated from young leaves following CTAB (Cetyl Trimethyl 
Ammonium Bromide) method as modified by Saghai-Maroof et al. 
[31] and analyzed by using SSR markers [32]— [33]. Markers 
spanning the whole genome at approximately 10 cM or closer 
intervals were used to assess parental polymorphism among the 
donor and the recipient line using 2.5% agarose gel. Out of 200 
SSR markers screened, 100 markers were polymorphic. To 
precisely estimate the size of alien chromosome segments, 
additional markers were screened in the regions harboring these 
introgressions. These polymorphic markers were used for analyz- 
ing the introgression and generating graphical genotypes of each 
introgression line using the software GGT 2.0 [34]. 

Results 

Agronomic performance of introgression line hybrids 
(ILHs) 

Analysis of variance revealed that mean squares attributable to 
genotypes were highly significant for all the traits, indicating 
considerable amount of genetic variation in the material used in 
this study (Table SI, Figure SI a-f). Out of the 318 test hybrids 
evaluated, some recorded a significant increase in yield compo- 
nents (Figure SI, a-f) while others, especially ones that had 
introgressions from O. glaberrima showed 100% spikelet sterility in 
some cases. Based on better yielding ability and phenotypic 
performance over the check hybrids, 17 hybrids were identified. 
These hybrids and their respective check hybrids had similar 
genetic background except for introgressed segments. Mean grain 
yield and percent standard heterosis of superior introgression line 
hybrids over their respective check hybrids are presented in 
Table 1. Four hybrids, viz., ILH299, ILH326, ILH411 and 
ILH435, were having Pusa44 with introgression enriched genome 
as male parent showed significantly higher yield over check hybrid 
PMS17A/Pusa44. This increase in standard heterosis ranged from 



37% to 40% over the check hybrid. ILH299 and ILH326 carried 
introgressions from O. rufipogon, whereas ILH41 1 and ILH435 had 
introgressions from O. nivara. In the second set of 13 hybrids, 
having enriched genome of PR 114, 10 had significant positive 
heterosis (22.6 to 44.3%) for grain yield over check hybrid 
PMS17A/PR114. ILH867 (PMS 1 7A/IL867) having introgres- 
sions from O. rufipogon (Acc. IRGC 104433) showed highest 
heterosis of 44.32%. These 17 hybrids were also superior over 
their respective recurrent parents, viz., Pusa44 and PR114, for 
yield and yield-contributing traits under study (data not shown). 
From this initial evaluation, it was clear that introgression 
segments from wild species in hybrid background contributed 
positively for yield improvement. 

In 2009, the 17 hybrids were again evaluated, along with their 
parents, including introgression lines and PMS17B (maintainer of 
PMS 17 A), two check hybrids, viz., PMS17A/Pusa44 and 
PMS17A/PR114, and recurrent parents (Pusa44 and PR114) for 
yield and yield-contributing traits. Analysis of variance revealed 
significant genotypic differences for all the traits except productive 
tillers per plant (Table SI). To assess the effect of alien 
introgression segments in homozygous state, the introgression 
lines were compared with their recurrent parents, i.e., Pusa44 and 
PR 114, for important yield-contributing traits (Table 2). Intro- 
gression lines, viz., IL299, IL326, IL41 1 and IL435, showed lower 
yield, whereas IL326 showed significantly fewer spikelets per 
panicle and grains per panicle than their recurrent parent Pusa44. 
Out of 13 introgression lines in PR1 14 genetic background, seven 
lines had marginal but non-significant increase in yield over 
PR114. IL867, IL894, IL921 and IL1459 showed significant 
increase for spikelets per panicle over PR1 14; however, they failed 
to show more grains per panicle than its recurrent parent because 
of lower spikelet fertility. None of introgression lines showed any 
improvement in grains per panicle and 1000-grain weight over 
PR 114. Results indicated that most of these introgression lines 
failed to surpass their recurrent parents for important yield- 
contributing traits. 

To assess the effect of introgression in heterozygous condition, 
the 17 ILs were crossed with CMS line PMS17A. The extent of 
heterosis in these hybrids for important agronomic traits over the 
parents, check hybrids and recurrent parents of introgressions lines 
is presented in Table 3. Heterosis for grain yield is attributable to 
simultaneous manifestation of heterosis for various yield compo- 
nents. Out of the four hybrids, having Pusa44-derived introgres- 
sion lines as one of the parent, ILH299 and ILH326 recorded 
significant positive heterosis to the tune of 30.4% and 74.3%, 
respectively, over their respective parental introgression lines 
IL299 and IL326. Both these hybrids recorded significant standard 
heterosis (>30%) over the check hybrid PMS 1 7A/Pusa44. 
Remaining two hybrids, viz., ILH41 1 and ILH435, failed to show 
any improvement over parents and checks for grain yield. Out of 
1 3 hybrids, having PR 1 1 4-derived introgression lines as one of the 
parent, five, viz., ILH867, ILH901, ILH921, ILH951 and 
ILH 1011, recorded significant positive heterosis over check hybrid 
PMS17A/PR114. ILH951 (PMS1 7A/IL951) recorded highest 
magnitude of heterosis for grain yield over IL951 (51.6%), 
PMS17A/PR114 (43.0%), PMS17B (18.4%) and PR114 
(22.8%). Except for ILH101 1, the remaining four hybrids showed 
a significant increase in grain yield over PR 114. No improvement 
in yield was observed because of low grains per panicle and low 
spikelet fertility in rest of eight hybrids. Five hybrids, viz., ILH299, 
ILH326, ILH867, ILH921 and ILH951, which gave higher grain 
yield than introgression line parent and check hybrids, were also 
seen to possess higher grain weight. ILH299, ILH326, ILH86 and 
ILH901 carried introgressions from O. rufipogon, whereas ILH921 
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Table 1. Mean grain yield and standard heterosis of 17 ILHs evaluated in the year 2008. 



Percent heterosis over respective 



Hybrid ID 


Cross 


Grain yield (t/ha) 


check hybrids 


ILH299 


PMS17A/IL299 [Pusa44/0. rufipogon IRGC100219//2*Pusa44] 


7.691 


37.1** 


ILH326 


PMS17A/IL326 [Pusa44/0. rufipogon IRGC100219//2*Pusa44] 


7.595 


35.4** 


ILH411 


PMS17A/IL411 [Pusa44/0. nivara (CR1001 1 3A//2*Pusa44] 


7.637 


36.2** 


ILH435 


PMS17A/IL435 [Pusa44/0. nivara (CR100142A//2*Pusa44] 


7.829 


39.6** 


ILH867 


PMS17A/IL867 [PR114/0. rufipogon IRGC104433//2*PR1 14] 


8.906 


44.3** 


ILH873 


PMS17A/IL873 [PR114/0. rufipogon IRGC104433//2*PR1 14] 


7.254 


26.2* 


ILH984 


PMS17A/IL894 [PR114/0. rufipogon IRGC104433//2*PR1 14] 


7.375 


19.5 


ILH896 


PMS17A/IL896 [PR114/0. rufipogon IRGC104433//2*PR1 14] 


7.702 


24.8* 


ILH900 


PMS17A/IL900 [PR114/0. rufipogon IRGC104433//2*PR1 14] 


7.568 


22.6* 


ILH901 


PMS17A/IL901 [PR114/0. rufipogon IRGC104433//2*PR1 14] 


8.102 


32.6** 


ILH921 


PMS17A/IL921 [PR114/0. nivara IRGC100593//2*PR114] 


8.097 


31.2** 


ILH924 


PMS17A/IL924 [PR114/0. nivara IRGC105695//2*PR114] 


7.377 


19.5 


ILH929 


PMS17A/IL929 [PR114/0. nivara IRGC100114//2*PR114] 


7.977 


29.2** 


ILH951 


PMS17A/IL951 [PR114/0. nivara (CR100142A//2*PR1 14] 


7.810 


19.1 


ILH1011 


PMS17A/IL101 1 [PR1 14/0. nivara (CR100140//2*PR1 14] 


8.391 


35.9** 


ILH1459 


PMS17A/IL1459 [PR114/0. glumaepatula (IR104387//2*PR1 14] 


8.034 


36.96** 


ILH1433 


PMS17A/IL1433 [PR114/0. longistaminata (IR1 04301 B//2*PR1 14] 


7.558 


23.35* 


Check Hybrid 1 


PMS17A/Pusa44 


5.627 




Check Hybrid 2 


PMS17A/PR114 


6.127 





*, **significant at P 
doi:10.1371/journal 



<0.05 and P<0.01, respectively. IRGC and CR refer to International Rice Germplasm Centre and CRRI Cuttack, accession numbers, respectively. 
pone.0096939.t001 



and ILH951 had introgressions from 0. nivara. Yield increase in 
these six hybrids was not associated with delayed flowering or 
maturity. In fact, almost all hybrids were earlier in flowering than 
the check hybrids, introgression lines and the recurrent parents by 
3 to 5 days (Table S2). Two hybrids (ILH867 and ILH921) were 
taller than check hybrids and PR114 but shorter than respective 
parental introgression lines IL867 and IL921. Remaining hybrids 
did not show any change in plant height. Thus, introgression for 
yield was not associated with change in plant architecture and 
prolonged crop duration in this study (Table S2). Genetic variation 
for productive tillers per plant was non-significant. However, for 
other yield-contributing traits, viz., panicle length and spikelets per 
panicle, ILH299, ILH326, ILH867, ILH873, ILH901, ILH921 
and ILH951 showed significant improvement over introgressions 
lines, check hybrids and respective recurrent parents. For panicle 
length, ILH921 (PMS17A/IL921) recorded 23% heterosis over 
male parent, whereas ILH951 (PMS17A/IL951) recorded 43.3% 
and 54.5% heterosis over IL951 and PR1 14, respectively. ILH921 
had highest panicle length among all the hybrids (Fig. 3). For 
grains per panicle, four hybrids, viz., ILH299, ILH326, ILH921 
and ILH951, recorded significant positive heterosis (>30%) over 
their respective introgression lines, check hybrid PMS17A/Pusa44 
(20-35%) and Pusa44 (24-38%). In 2008, out of 318 hybrids 
tested, 60 hybrids having introgressions from O. glaberrima showed 
almost 100% spikelet sterility. Similarly, in 2009, almost all 
hybrids had lower spikelet fertility than recurrent parents except 
for hybrid ILH921, which expressed >10% significant heterosis 
over IL921, PMS17B and PMS17A/PR1 14. Spikelet sterility was 
inversely related to yield in hybrids. These hybrids were also tested 
for pollen fertility. High yielding hybrids, viz., ILH299, ILH326, 
ILH921 and ILH951, had higher pollen fertility than the check 
hybrids. From these results, it can be seen that introgressed 



segments in homozygous state failed to show improvement in grain 
yield in the introgression lines. 

Marker analysis and generation of graphical genotypes 

Introgression lines of six promising hybrids, viz., IL299, IL326, 
IL867, IL901, IL921 and IL951, which exhibited significant 
heterosis over their respective check hybrids, were used for the 
identification of alien segments responsible for enhanced perfor- 
mance of yield-related traits. A total of 200 SSR markers spanning 
all 12 chromosomes were used for characterizing each introgres- 
sion line. Introgression lines, viz., IL299, IL326, IL867 and IL901, 
were in BC 2 F 6 generation, whereas IL921 and IL951 were in 
BC 2 F 8 generation. Introgression lines IL299 and IL326, which 
shared a common parentage, had 16 and 10 marker introgres- 
sions, (7.7% and 5.2%, respectively) from 0. rufipogon. Introgres- 
sion lines IL867 and IL901, which also shared a common 
parentage, had 3.4% and 2.5% genome introgressions from O. 
rufipogon. In introgression line IL921, only five markers (2.2%) 
showed introgressions from O. nivara, whereas line IL951 showed 
12 marker introgressions (6.7%) from 0. nivara. Details of 
introgressions in six introgression lines are presented in Table 
S3. Depending upon the introgression pattern of six introgression 
lines, graphical genotypes for each line were generated (Fig. 4). 

Associating yield and yield components with alien 
introgressed regions 

Yield increase in hybrids, viz., ILH299, ILH326, ILH867, 
ILH901, ILH921 and ILH951, was caused by increase in panicle 
length, spikelets per panicle, grains per panicle and 1000-grain 
weight. Hybrids ILH299, ILH326, ILH867 and ILH901 showed 
increased 1000-grain weight over their respective introgression 
lines IL299, IL326, IL867 and IL901. Comparison of IL299, 
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IL326, IL867 and IL901 revealed common introgressed regions 
on chromosomes 3 and 1 1 . All the ILs had common introgressions 
around 191-202 cM region on chromosome 3 and 40-102 cM 
region on chromosome 11 (Fig. 4). Likewise, IL921 and IL951 
derived from crosses with 0. nivara had common introgressions 
around 31 cM region on chromosome 5. Introgression line IL299 
had 0. rufipogon allele at RM514, RM148 and RM442, whereas 
IL326 had RM416 on chromosome 3. IL867 carried 0. rufipogon 
aUeles at RM448, RM514, RM570, RM148 and RM1 14, whereas 
IL901 also had introgression for marker RM448 (Fig. 4). In the 
present study, ILH95 1 showed higher magnitude of heterosis for 
1000-grain weight than IL951 and the checks. Introgression line 
IL951 had 0. nivara alleles at RM222 (11.3 cM) on chromosome 
10 and RM254 (1 10.0 cM) on chromosome 1 1 (Fig. 4; Fig. S3). O. 
rufipogon and O. nivara alleles at these marker loci could be 
associated with increased 1000-grain weight in introgression line 
hybrids. ILH299 and ILH326 showed improvement in spikelets 
per panicle, grains per panicle and spikelet fertility over IL299 and 
IL326, respectively. Both the introgression lines had common 
introgression at marker loci RM508 (0.0 cM) on chromosome 6, 
RM407 (5.7 cM) on chromosome 8, RM202 (54.0 cM) to RM21 
(85.0 cM) on chromosome 11 and RM19 (20.9 cM) on chromo- 
some 12 (Fig. 4; Fig. S2). These introgression lines showed 
introgression of O. rufipogon chromatin at RM84, RM220 and 
RM47 3 A on chromosome 1 and at RM2 1 on chromosome 1 1 
(Fig. 4). ILH921 recorded higher yield than IL921, PMS17A/ 
PR 1 1 4 and parents because of increased panicle length, spikelet 
fertility, spikelets per panicle and grains per panicle. Molecular 
characterization of IL921 revealed introgressions of O. nivara 
chromatin at RM302 (147.8 cM) and RM472 (171.6 cM) on 
chromosome 1 and RM307 (0.0 cM) and RM401 (8.5 cM) on 
chromosome 4. ILH951 was the highest yielding among all 
hybrids because of highest spikelets per panicle, grains per panicle, 
spikelet fertility and 1000-grain weight. Parental introgression line 
(IL95 1 ) of this hybrid had introgression from O. nivara in seven 
chromosomes. Introgression was observed for RM109 (4.8 cM) on 
chromosome 1, RM289 (56.7 cM) on chromosome 5, RM584 
(20.2 cM) and RM 351 (73.3 cM) to RM505 (78.6 cM) on 
chromosome 7, RM38 (28.0 cM) on chromosome 8, RM474 
(0.0 cM) to RM222 (11.5 cM) on chromosome 10, and RM206 
(102.9 cM) to RM254 (1 10.0 cM) on chromosome 1 1 (Fig. 4). The 
observed improvement in yield-contributing traits in this hybrid 
could be attributable to introgressions from O. nivara. In addition 
to yield components, days to 50% flowering and plant height 
decides the suitability of a hybrid for production. Introgression 
lines IL867 and IL901 were earlier than recurrent parent PR1 14, 
but their hybrids ILH867 and ILH901 took more days to mature 
than parental introgression lines. O. rufipogon alleles at marker locus 
RM287 on chromosome 11 in IL867 and IL901 might show 
increased effect on days to 50% flowering in their hybrids. 
Introgression line IL921 showed a significant increase in plant 
height over recurrent parent PR114, but its hybrid ILH921 had 
reduced plant height over IL921. O. nivara alleles at RM472 on 
chromosome 1 and RM307 on chromosome 4 could be 
responsible for increased plant height in IL921, but when these 
alleles are present in heterozygous condition in ILH92 1 , there was 
reduction in plant height. Similarly, ILH95 1 showed a significant 
reduction in days to 50% flowering and days to maturity relative to 
IL951 and PMS17A/PR114. 

Discussion 

Though rice has a wide genetic base, the germplasm with which 
breeders are working, in general, has limited variability, which 
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Figure 3. Field photograph of IL921 (left), ILH921 (middle) and 

doi:10.1371/journal.pone.0096939.g003 

probably is the major bottleneck in improving productivity. 
Because breeders, in general, are not currently using unadapted 
germplasm in their breeding programmes, pre-breeding must be 
strengthened to widen the genetic base of the working germplasm. 
Incorporation of 'hidden genes' in improved rice germplasm, 
through development and utilization of alien introgression lines for 
productivity enhancement, has emerged as a good breeding 
strategy [4]. The common wild rice species, 0. rufipogon and 0. 
nivara, are known as the ancestors of cultivated rice and have been 
proven to be a valuable gene pool for rice genetic improvement 
[35] . Sun et al. [36] studied the genetic diversity of 0. rufipogon and 



right hand side showing panicle characteristics. 

0. sativa from more than 10 Asian countries and found that the 
genetic diversity of 0. rufipogon was much higher than that of 
cultivated rice. A large number of alleles present in O. rufipogon 
could not be found in cultivated rice. Utilizing these hidden genes 
for the genetic enrichment of rice cultivars/hybrids could bring a 
revolutionary change in productivity improvement. 

Introgression lines (ILs) serve as a valuable tool for genetic 
studies and genome analysis. A common accession of 0. rufipogon 
(IRGC 105491) was used in the development of alien ILs in most of 
the previous studies [14]-[16], [19], [37]-[39]. In the present 
study, out of 17 ILs used for generating Fj hybrids, six were 
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Figure 4. Graphical genotypes of six alien introgression lines generated after analyzing these with polymorphic SSR markers. 

Regions in blue are homozygous alien segments and gray are heterozygous alien segments introgressed from O. rufipogon in IL299, IL326, IL867 and 
IL901 and O. nivara in IL901 and IL951 . Numbers on the right of linkage group are the cM distances as per Temnykh et al [32], Numbers at the bottom 
are the chromosomes. Details of the ILHs in parentheses are presented in Table 1. 
doi:1 0.1 371 /journal. pone.0096939.g004 
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Table 4. Correspondence of reported introgressions from 0. rufipogon and 0. nivara with the already reported quantitative trait 
locus/loci. 





Putative 
traits 


ILs 


Donor species 


Introgression 
reported in 
chromosome 


Marker associated 
with introgressed 
regions 


Marker associated with QTL(s) 
and the species as reported in 
previous studies 


Reference 


Panicle length 
(cm) 


IL921 


0. nivara 


1 


RM472 


RM472 (p//) O. rufipogon, RM315 (pli.l) 
O. rufipogon, RM315 {pll.i) O. rufipogon 


[40,16,37] 


Spikelets per 
panicle 


IL299 


O. rufipogon 


1 


RM220 


RG532 (sppi.2) 0. rufipogon, RM246 (sppi.2) 
O. rufipogon,RM220 [spp1) O. glumaepatula 


[14,40,41] 




IL326 


O. rufipogon 


1 


RM473A 


RM237 {sppl.2) O. sativa cv. Jefferson 


[16] 




IL921 


O. nivara 


4, 1, 5 


RM307, RM302, RM413 


RM307 (hsp4) O. rufipogon, RZ730 (sppl.1) 
O. rufipogon, RM194 (snp5.1) O. rufipogon 


[20,15,42] 


Grains per 
panicle 


IL299 


O. rufipogon 


1 


RM220, RM84 


RM220 (gpll.l) O. rufipogon, RM220 
(ggpi) 0. glumaepatula, RM84 iqGPLI) 
O. rufipogon 


[37,41,43] 




IL951 


O. nivara 


7, 10 


RM11, RM222 


RM11 (qFG7) O. rufipogon, RM222 (hgpW) 
O. rufipogon, C153A-RM222 (gpW) 
O. sotivo 


[43,20,44] 




IL921 


O. nivara 


5 


RM413 


RM194 (gnS.I) O. rufipogon 


[42] 


Spikelet 
fertility (%) 


IL921 


O. nivara 


1, 4 


RM302, RM472, RM307 


RM212 (pss/.!) 0. rufipogon, RM315 (sfl.i) 
O. rufipogon, RM315 [pssl.i) O. sativa cv. 
Jefferson, RM307 {hspp4) 0. rufipogon 


[37,42,16,20] 




IL951 


O. nivara 


7, 5 


RM11, RM289 


RM1 1 (qss7) O. rufipogon, RM289 (hsspS) 
0. rufipogon 


[43,20] 


1000 grain 
weight (g) 


IL299 


O. rufipogon 


3, 11 


RM442, RM21 


RM227 (gw3.1) 0. rufipogon, RM21 (hgwll) 
0. rufipogon, RG2-RM21 (gwUb) 0. sativa 


[15,20,44] 




IL323 


O. rufipogon 


11 


RM21 


RM21 (hgwl 1) 0. rufipogon 


[20] 




IL867 


O. rufipogon 


3 


RM514, RM570, RM180, 
RM114 


RM227 (gw3.1) 0. rufipogon 


[16] 




IL951 


O. nivara 


10, 11 


RM222, RM254 


RM222 {hgwIO) 0. rufipogon, 
C153A-RM222 {gwlO) 0. sativa, 
G89-G1084 (gwlO) O. sativa, RM254 
{gwii.2) O. rufipogon 


[20,44,45,15] 


Grain Yield 
(t/ha) 


IL299 


O. rufipogon 


1 


RM220 


RM220 (y/dJ.I) 0. rufipogon, RM220 (gypl) 
0. glumaepatula, RM220 0. glumaepatula 


[37,41,23] 




IL951 


O. nivara 


8, 2 


RM38, RM109 


RM38 (y/dS.3) 0. rufipogon, RM236 (hyp2a) 
0. rufipogon 


[42,20] 



Markers in bold are common regions associated with yield and yield components in various studies. 
doi:1 0.1 371 /joumal.pone.0096939.t004 



derived by crossing two indica rice cultivars (Pusa44 and PR 11 4) 
with two different accessions each of 0. rufipogon (IRGC100219 
and IRGC 104433) and 0. nivara (IRGC 100593 and CR100142A). 
These were used for identification of alien segments responsible for 
enhanced magnitude of heterosis for yield and its component 
traits. In these ILs, 37 markers confirmed the presence of 32 alien 
segments dispersed across the genome. Sixteen alien segments, 
associated with 19 markers, have been located in the same or 
similar regions where QTL for yield and its component traits were 
mapped in previous studies [14]-[16], [20], [23], [37], [40]-[45] 
(Table 4). ILHs developed from these alien ILs, carrying 
introgressions from 'A' genome donor species, demonstrated the 
scope of utilization of alien introgressions in heterosis breeding. 
Results from two years of evaluation revealed that hybrids having 
introgressions from 0. rufipogon (ILH299, ILH326, ILH867and 
ILH901) and 0. nivara (ILH921 and ILH951) in their genome, 
showed significant heterosis for grain yield, panicle length, 
spikelets per panicle, grains per panicle and 1000-grain weight 
over introgression lines, check hybrids and hybrids generated by 
crossing the cms line and the recurrent parent. Superiority of 
introgression lines over recurrent parents have been reported by 
several workers, but the studies on heterosis involving alien 



introgression lines are quite limited. Xiao et al. [14] identified 
QTL for panicle length (pll.l, pll.2, pl2.1, pl4.1, pW.l, pl9.1& 
pl9.1), spikelets per panicle (sppl.l, sppl.2 & spp9.1), grains per 
panicle (gppl .1, gpp8.2 &gppl2.i), percent seed set {pss2.1& pss4.1), 
and 1000-grain weight (gw4.1, gwS.l, gw9.1, gwll.l & gwl2.1) 
from 0. rufipogon in interspecific BC 2 testcross population. O. 
rufipogon allele at two QTL, namely, jldl.l and yld2.1, was 
associated with 18 and 17 per cent grain yield increment, 
respectively. These two QTL were later transferred to the elite 
restorer lines Minghui 63, 93-11 and Ce64-7 through marker 
assisted selection and their functions were confirmed by significant 
yield improvement in rice hybrids [45] . In another study, involving 
F! hybrids between a cross of 0. rufipogon-devived ILs and recurrent 
parent Guichao2, 30 heterotic loci associated with six yield-related 
traits were reported by Luo et al. [20] . In tomato, Gur and Zamir 
[46] successfully increased fruit yield using ILs carrying a pyramid 
of three independent yield-promoting genomic regions introduced 
from wild species Solanum pennellii. Yield of hybrids developed using 
these pyramided ILs was more than 50% higher than that of a 
leading commercial variety. This study demonstrated the merit of 
'exotic alleles' in improving crop productivity. 



PLOS ONE | www.plosone.org 



9 



June 2014 | Volume 9 | Issue 6 | e96939 



Yield QTL from Wild Rice 



Most of the markers associated with the alien introgressions, in 
this study, are in correspondence with the already reported QTL 
for yield-contributing traits from 0. rufipogon [14]— [16], [37], [40], 
[42]— [43] (Table 4). In accordance with the earlier reports, results 
of the present study also showed the positive association of alien 
segments with the yield-contributing traits and their utilization in 
hybrid development. Two experimental hybrids, viz., ILH921 
(PMS17A/IL921) and ILH951 (PMS17A/IL951), developed by 
using PR1 14-derived ILs, carried 0. nivara chromosomal segment 
in their genome. ILH921 showed increased yield over parent 
IL921 and check hybrid through yield-contributing traits. IL921 
had 0. nivara introgressions around markers RM302, RM472 in 
chromosome 1 and around marker RM307 in chromosome 4, 
which are same as where 0. rufipogon alleles have been observed to 
confer positive effect on panicle length, spikelets per panicle and 
spikelet fertility in previous studies [14], [16], [20], [37], [40], [42] 
(Table 4). ILH951 was highest yielding among all the hybrids. 
Pollen parent of this hybrid, IL95 1 , had introgressions from 0. 
nivara in nine regions on seven chromosomes. Marker loci RM1 1, 
RM222, RM254, RM289 and RM109 on chromosomes 7, 10, 1 1, 
5 and 2, respectively, exhibited presence of O. nivara chromatin. 
The same markers were reported to be linked to QTL transferred 
from O. rufipogon for grains per panicle, spikelet fertility, 1000-grain 
weight and yield per plant [15], [20], [42], [43] (Table 4). Markers 
responsible for improvement in yield traits through incorporation 
of genomic regions from 0. nivara in the present study were also 
reported to be linked with QTL derived from O. rufipogon. Though 
these two species have distinct phenotypes because of their 
different life cycles and breeding behavior in their natural habitats, 
because of their evolutionary path they share some common/ 
conserved genomic regions for yield component traits, as also 
observed in this study. O. rufipogon alleles seem to be effective across 
recipient genotypes and evaluation sites and O. nivara alleles 
showed synteny with the alleles of O. rufipogon. This synteny 
between genomic regions derived from both the alien species is 
well demonstrated by SSRs. Occurrence of intermediate popula- 
tions between O. nivara and O. rufipogon types in some habitats and 
introgressions among these wild species and between cultivated 
and wild species [47]-[49] may be the reason for such homology. 
The complete information on the extent of introgression between 
these two wild species is still an enigma [50]-[53]. 

Pollen fertility of hybrids is important for realization of superior 
yields. Higher pollen fertility (>10%) in ILH299, ILH326, 
ILH901, ILH921 and ILH951 than that of check hybrids indicates 
the role of alien species-derived modifier genes present in ILs. 

Introgression lines, viz., IL299, IL326, IL867, IL921 and IL951, 
did not show any improvement in grains per panicle, 1000-grain 
weight and grain yield when compared with recurrent parents, but 
their hybrids showed significant improvement over check hybrids 
and recurrent parents. Alien segments in heterozygous state 
showed an increasing effect for these traits (overdominance). It 
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seems that in heterozygous condition, the recipient complement is 
conserved and the minor negative effects associated with alien 
segments are masked. Introgression in this set of ILs has been 
recorded in several chromosomal regions and some of these are 
large in size. These QTL, thus, need to be mapped for their 
precise transfer into other elite parental lines, as the superiority of 
the alien segment in heterozygous state can be readily exploited in 
hybrid rice breeding. The strength of utilization of related wild 
species in improving level of heterosis is well demonstrated 
through this experiment. 

Supporting Information 

Figure SI (a-f): Frequency distributions of 318 ILHs 
for yield and yield contributing traits. 

(EPS) 

Figure S2 Ethidium bromide stained DNA amplifica- 
tion profile of parents and introgression lines using 
microsatellite markers. Lanes are: M=100bp ladder, 

1 =Pusa44, 2 = 0. rufipogon (Acc. IR 104852), 3 and 4 = introgres- 
sion lines (IL299 and IL326), C = negative control. 

(EPS) 

Figure S3 Ethidium bromide stained DNA amplifica- 
tion profile of parents and introgression lines using 
microsatellite markers. Ethidium bromide stained DNA 
amplification profile of parents and introgression lines using 
microsatellite markers. Lanes are: M = 100 bp ladder, 1 = PR1 14, 

2 = 0. nivara (Acc. IR100142A), 3 = introgression line (IL951), 
C = negative control. 

(EPS) 

Table SI Mean square values for different agronomic 
traits based on ANOVA of field trials of introgression 
line hybrids. 

(DOC) 

Table S2 Mean values for some important agronomic 
traits and extent of heterosis (%) over parents, check 
hybrids and commercial checks in 17 hybrids. 

(DOC) 

Table S3 Proportion of microsatellite markers that 
observed introgression from O. rufipogon (IL299) and 
O. nivara (IL951). 

(DOC) 

Author Contributions 

Conceived and designed the experiments: KS KBG. Performed the 
experiments: KBG NS DB RK KS. Analyzed the data: KBG NS NSB 
TSB DB KS. Contributed reagents/materials/analysis tools: KS TSB NSB 
NS. Wrote the paper: KBG KS NS NSB TSB. Helped in development of 
introgression lines: DB KS. Guided the study: NSB TSB. 



8. Khush GS, Bacalangeo E, Ogawa T (1990) A new gene for resistance to 
bacterial blight from 0. bngis/aminala. Rice Genetics Newsletter 7: 121-122. 

9. Zhang Q, Zhao BY, Zhao KJ, Wang CL, Yang WC, ct al. (2000) Identifying 
and mapping a new gene Xa23(t) for resistance to bacterial blight {Xanthornonas 
oryzae pv. oryzae) from 0. rufipogon. Acta Agronomica Sinica 26: 536-542. 

10. Bhasin H, Bhatia D, Raghuvanshi S, LoreJS, Sahi GK, et al. (2012) New PCR- 
based sequence-tagged site marker for bacterial blight resistance gene Xa38 of 
rice. Molecular Breeding 30: 607-611. 

11. Frcy KJ, Cox TS, Rodgcrs DM, Cox PB (1984) Increasing cereal yields with 
genes from wild and weedy species. Genetics: New Frontiers Proc XV Intern 
Cong Genet Oxford and IBH Publishing Co New Delhi 4: 51-68. 

12. Tanksley SD, Nelson JC (1 996) Advanced backcross QTL analysis: a method for 
the simultaneous discovery and transfer of valuable QTLs from unadapted 



PLOS ONE | www.plosone.org 



10 



June 2014 | Volume 9 | Issue 6 | e96939 



Yield QTL from Wild Rice 



germplasm into elite breeding lines. Theoretical and Applied Genetics 92: 191- 
203. 

1 3. Xiao J, Grandillo S, Ahn SN, McCouch S, Tanksley SD (1996) Genes from wild 
rice improves yield. Nature 384: 223-224. 

14. Xiao J, Li J, Grandillo S, Yuan L, Tanksley SD, et al.(1998) Identification of 
traits improving quantitative trait loci alleles from wild rice relative Oryza 
rufipogon. Genetics 150: 899-909. 

15. Moncada P, Martinez CP, Borrero J, Chatel M, Gauch JH, et al. (2001) 
Quantitative trait loci for yield and yield components in an Oryza sativaxOryza 
rufipogon BC 2 p2 population evaluated in an upland environment. Theoretical and 
Applied Genetics 102: 41-52. 

16. Thomson MJ, Tai TH, McClung AM, Lai XH, Hinga ME, et al. (2003) 
Mapping quantitative trait loci for yield, yield components and morphological 
traits in an advanced baekcross population between Oryza rufipogon and the Oryza 
saliva eultivar Jefferson. Theoretical and Applied Genetics 107: 479—493. 

17. He G, Luo X, Tian F, Li K, Zhu Z, et al. (2006) Haplotypc variation in structure 
and expression of a gene cluster associated with a quantitative trait locus for 
improved yield in rice. Genome research 16: 618-626. 

18. Xie X, Jin F, Song MH, SuhJP, Hwang HG, et al. (2008) Fine mapping of a 
yield- enhancing QTL cluster associated with transgressive variation in an Oryza. 
sativaxO. rufipogon cross. Theoretical and Applied Genetics 116: 613-622. 

19. Jin FX, Kim DM, Ju HG, Ahn SN (2009) Mapping quantitative trait loci for 
awnness and yield component traits in isogenic lines derived from an Oiyza 
saliva/ 0. rufipogon cross. Journal of Crop Science and Biotechnology 12: 9—15. 

20. Luo X, Wu S, Tian F, Xin X, Zha X, et al. (201 1) Identification of heterotic loci 
associated with yield-related traits in Chinese common wild rice (Oryza rufipogon 
Griff.). Plant Science 181: 14-22. 

21. Aluko G, Martinez C, Tohme J, Castano C, Bergman C, et al. (2004) QTL 
mapping of grain quality traits from the interspecific cross Oryza sativaxO. 
glaberrima. Theoretical and Applied Genetics 109: 630-639. 

22. Yoon DB, Kang KH, Kim HJ, Ju HG, Kwon SJ, et al. (2006) Mapping 
quantitative trait loci for yield components and morphological traits in an 
advanced baekcross population between Oryz.a grandiglumis and the 0. sativa 
japonica eultivar Hwaseongbyeo. Theoretical and Applied Genetics 1 12: 1052— 
1062. 

23. Rangel PN, Brondani RP, Rangel PH, Brondani C (2008) Agronomic and 
molecular characterization of introgression lines from the interspecific cross 
Oryza sativa (BG90-2) x Oryza glumaepatula (RS-16). Genet Mol Res 7: 184-195. 

24. Chen Z, Hu F, Xu P, Li J, Deng X, et al. (2009) QTL analysis for hybrid sterility 
and plant height in interspecific populations derived from a wild rice relative, 
Oryza longislaminala. Breeding science 59: 441—445. 

25. Linh LH, Hang NT, Jin FX, Kang KH, Lee YT, et al. (2008) Introgression of a 
quantitative trait locus for spikelets per panicle from Oryza minuta to the 0. sativa 
eultivar Hwaseongbyeo. Plant Breeding 127: 262—267. 

26. Eshed Y, Zamir D (1995) An introgression line population of lycopersicon pennellii 
in the cultivated tomato enables the identification and fine mapping of yield- 
associated QTL. Genetics 141: 1147-1162. 

27. Monforte AJ, Tanksley SD (2000) Fine mapping of a quantitative trait locus 
(QTL) from Lycopersicon hirsutum chromosome 1 affecting fruit characteristics and 
agronomic traits: breaking linkage among QTLs affecting different traits and 
dissection of heterosis for yield. Theoretical and Applied Genetics 100: 471—479. 

28. Semel Y, Nissenbaum j, Menda N, Zinder M, Krieger U, et al. (2006) 
Overdominant quantitative trait loci for yield and fitness in tomato. Proc Nad 
Acad Sei USA 103: 12981-12986. 

29. Standard Evaluation System for Rice (2002) International Rice Research 
Institute, Manila, Philippines. 

30. SPAR 2.0: Statistical Package for Agricultural Research Ver. 2.0, Indian 
Agricultural Statistics Research Institute, New Delhi, India, (http://iasri.res.in/ 
spar) 

3 1 . Saghai-Maroof MA, Soliman KM, Jorgensen AR, Allard RW ( 1 984) Ribosomal 
DNA spacer length polymorphism in barley: Mendelian inheritance, chromo- 
somal location and population dynamics. Proc Natl Acad Sei USA 81: 8014- 
8018. 

32. Temnykh S, Park WD, Ayres N, Cartinhour S, Hauck N, et al. (2000) Mapping 
and genome organization of microsatcllitc sequences in rice (Oryza sativa L.). 
Theoretical and Applied Genetics 100: 697-712. 



33. McCouch SR, Teytclman L, Xu Y, Lobos KB, Clare K, et al. (2002) 
Development and mapping of 2240 new SSR markers for rice (Oryza sativa L.). 
DNA research 9: 199-207. 

34. Van Berloo R (2008) GGT 2.0: versatile software for visualization and analysis of 
genetic data. Journal of Heredity 99: 232-236. 

35. Khush GS (1997) Origin, dispersal, cultivation and variation of rice. Plant 
Molecular Biology 35: 25-34. 

36. Sun CQ, Wang XK, Li ZC, Yoshimura A, Iwata N (2001) Comparison of the 
genetic diversity of common wild rice (Oryza rufipogon Griff.) and cultivated rice 
(0. sativa L.) using RFLP markers. Theoretical and Applied Genetics 102: 157- 
162. 

37. Septiningsih EM, Prasetiyono J, Lubis E, Tai TH, Tjubaryat T, et al. (2003) 
Identification of quantitative trait loci for yield and yield components in an 
advanced baekcross population derived from the Oryza sativa variety IR64 and 
the wild relative 0. rufipogon. Theoretical and Applied Genetics 107: 1419-1432. 

38. Cheema KK, Bains NS, Mangat GS, Das A, Vikal Y, et al. (2008) Development 
of high yielding IR64 xOnza rufipogon (Griff.) introgression lines and identifica- 
tion of introgrcsscd alien chromosome segments using SSR markers. Euphytica 
160: 401-409. 

39. Tan L, Liu F, Xue W, Wang G, Ye S, et al. (2007) Development of Oryza 
rufipogon and 0. sativa Introgression Lines and Assessment for Yield-related 
Quantitative Trait Loci. Journal of Integrative Plant Biology 49: 871-884. 

40. Lee SJ, Oh CS, SuhJP, McCouch SR, Ahn SN (2005) Identification of QTLs 
for domestication-related and agronomic traits in an Oryza sativaxO. rufipogon 
BC1F7 population. Plant Breeding 124: 209-219. 

41. Brondani C, Rangel P, Brondani R, Ferreira M (2002) QTL mapping and 
introgression of yield-relalcd traits from Onrji glumaepatula to cultivated rice 
(Oryza sativa) using microsatcllitc markers. Theoretical and Applied Genetics 104: 
1192-1203. 

42. Marri PR, Sarla N, Reddy LV, Siddiq EA (2005) Identification and mapping of 
yield and yield related QTLs from an Indian accession of Oryza rufipogon. BMC 
Genetics 6: 633-646. 

43. Tian F, Fu Q, Zhu ZF, Fu YC, Wang XK, et al. (2006) Construction of 
introgression lines carrying wild rice (Oryza rufipogon Griff.) segments in cultivated 
rice (Oryza sativa L.) background and characterization of introgressed segments 
associated with yield-related traits. Theoretical and Applied Genetics 1 12: 570— 
580. 

44. Hua J, Xing Y, Wu W, Xu C, Sun X, et al. (2003) Single-locus heterotic effects 
and dominance by dominance interactions can adequately explain the genetic 
basis of heterosis in an elite rice hybrid. Proc Natl Acad Sci USA 100: 2574— 
2579. 

45. Luo LJ, Li ZK, Mei HW, Shu QY, Tabien R, et al. (2001) Overdominant 
epistatic loci are the primary genetic basis of inbreeding depression and heterosis 
in rice. II. Grain yield components. Genetics 158: 1755-1771. 

46. Gur A, Zamir D (2004) Unused natural variation can lift yield barriers in plant 
breeding. PLoS Biology 2: 1610-1615. 

47. Oka HI (1988) Origin of cultivated rice: Japan Scientific Societies Press. 

48. Deng QY, Yuan LP, Liang FS, Li J, Wang LG, et al. (2004) Studies on yield- 
enhancing genes from wild rice and their marker-assisted selection in hybrid rice. 
Hybrid Rice 19: 6-10. 

49. Song Z, Zhu W, RongJ, Xu X, Chen J, et al. (2006) Evidences of introgression 
from cultivated rice to Oryza rufipogon (Poaccae) populations based on SSR 
fingerprinting: implications for wild rice differentiation and conservation. 
Evolutionary Ecology 20: 501-522. 

50. Vaughan DA, Morishima H, Kadowaki K (2003) Diversity in the Oryza genus. 
Current Opinion in Plant Biology 6: 139-146. 

51. Kuroda Y, Sato YI, Bounphanousay C, Kono Y, Tanaka K (2007) Genetic 
structure of three Oryza AA genome species (0. rufipogon, 0. nivara and 0. sativa) 
as assessed by SSR analysis on the Vientiane Plain of Laos. Conservation 
Genetics 8: 149-158. 

52. Sang T, Ge S (2007) Genetics and phylogenetics of rice domestication. Current 
Opinion in Genetics & Development 17: 533—538. 

53. Zhou HF, Zheng XM, Wei RX, Second G, Vaughan DA, et al. (2008) 
Contrasting population genetic structure and gene flow between Oryza rufipogon 
and Oryza nivara. Theoretical and Applied Genetics 117: 1181-11 89. 



PLOS ONE | www.plosone.org 



11 



June 2014 | Volume 9 | Issue 6 | e96939 



